Submitted to Physical Review Letters 



Two-scale structure of the electron dissipation region during collisionless magnetic 

reconnection 

M. A. Sha£| 

Department of Physics & Astronomy, 217 Sharp Lab, University of Delaware, Newark, DE 19716 



J. F. Drake, M. Swisdak 
University of Maryland, College Park, MD, 20742 
(Dated: February 5, 2008) 

Particle in cell (PIC) simulations of collisionless magnetic reconnection are presented that demon- 
strate that the electron dissipation region develops a distinct two-scale structure along the outflow 
direction. The length of the electron current layer is found to decrease with decreasing electron 
mass, approaching the ion inert ial length for a proton-electron plasma. A surprise, however, is that 
the electrons form a high- velocity outflow jet that remains decoupled from the magnetic field and 
extends large distances downstream from the x-line. The rate of reconnection remains fast in very 
large systems, independent of boundary conditions and the mass of electrons. 

PACS numbers: Valid PACS appear here 



Magnetic reconnection drives the release of magnetic 
energy in explosive events such as disruptions in labo- 
ratory experiments, magnetic substorms in the Earth's 
magnetosphere and flares in the solar corona. Recon- 
nection in these events is typically collisionless because 
reconnection electric fields exceed the Dreicer runaway 
field. Since magnetic field lines reconnect in a boundary 
layer, the "dissipation region" , whose structure may limit 
the rate of release of energy, understanding the structure 
of this boundary layer and its impact on reconnection 
is critical to understanding the observations. Because 
of their ability to carry large currents the dynamics of 
electrons continues to be a topic of interest. Early sim- 
ulations of reconnection suggested that the rate of re- 
connection was not sensitive to electron dynamics p], Q 
and this insensitivity was attributed to the coupling to 
whistler dynamics at the small spatial scales of the dis- 
sipation region 0, 0] • The results of more recent kinetic 
PIC simulations have called into question these results 
by suggesting that the electron current layer stretches 
along the outflow direction and the rate of reconnection 
drops The fast rates of reconnection obtained from 

earlier simulations N, 13, III were attributed to the influ- 
ence of periodicity[5[. 

We present particle-in-cell (PIC) simulations with var- 
ious electron masses and computational domain sizes and 
an analytic model that demonstrate that collisionless re- 
connection remains fast even in very large collisionless 
systems. The reconnection rate stabilizes before the pe- 
riodicity of the boundary conditions can impact the dy- 
namics. The electron current layer develops a distinct 
two-scale structure along the outflow direction that had 
not been identified in earlier simulations. The out-of- 
plane electron current driven by the reconnection elec- 
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trie field has a length that decreases with the electron 
mass, scaling as (m e /mj) 3 / 8 , which extrapolates to about 
an ion inertial length di — c/uo P i for the electron-proton 
mass ratio. The surprise is that a jet of outflowing elec- 
trons with velocity close to the electron Alfven speed CAe 
extends up to several 10's of di from the x-line. Remark- 
ably, the electrons are able to jet across the magnetic field 
over such enormous distances because momentum trans- 
port transverse to the jet effectively "blocks" the flow of 
the out-of-plane current in this region. The momentum 
transport causing this "current blocking" effect has the 
same source (the off diagonal pressure tensor but is 
much stronger than that which balances the reconnection 
electric field at the x-line. 

Our simulations are performed with the particle-in-cell 
code p3d [1, 0]. The results are presented in normal- 
ized units: the magnetic field to the asymptotic value of 
the reversed field, the density to the value at the cen- 
ter of the current sheet minus the uniform background 
density, velocities to the Alfven speed va, lengths to 
the ion inertial length di, times to the inverse ion cy- 
clotron frequency fi^ 1 , and temperatures to vfiiv\. We 
consider a system periodic in the x — y plane where flow 
into and away from the x-line are parallel to y and x, 
respectively. The reconnection electric field is parallel 
to z. The initial equilibrium consists of two Harris cur- 
rent sheets superimposed on a ambient population of uni- 
form density. The reconnection magnetic field is given by 
B x = tanh[(?/ - L y /A)/w ] - tanh[(?/ - 3L y /4)/w } - 1, 
where wo and L y are the half- width of the initial current 
sheets and the box size in the y direction. The electron 
and ion temperatures, T e = 1/12 and Ti = 5/12, are ini- 
tially uniform. The initial density profile is the usual Har- 
ris form plus a uniform background of 0.2. The simula- 
tions presented here are two-dimensional, i.e., d/dz — 0. 
Reconnection is initiated with a small initial magnetic 
perturbation that produces a single magnetic island on 
each current layer. 

We have explored the dependence of the rate of recon- 
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FIG. 1: (color online). Reconnection electric field versus time: 
(a) 204.8 x 102.4, (b) 102.4 x 51.2, (c) 51.2 x 25.6. w is the 
initial current sheet width. 

nection on the system size in a series of simulations with 
three different system sizes and three different mass ra- 
tios. For mi/m e = 25, the grid scale A = 0.05 and the 
speed of light c = 15. For rrii/m e = 100, A = 0.025 and 
c 20. For mi/m e = 400, A 0.0125 and c 40. The 
reconnection rate versus time is plotted for our simula- 
tions in Fig. [H The reconnection rate is determined by 
taking the time derivative of the total magnetic flux be- 
tween the x-line and the center of the magnetic island. 
The rate increases with time, undergoes a modest over- 
shoot that is more pronounced in the smaller domains, 
and approaches a quasi-steady rate of around 0.14, in- 
dependent of the domain size. Earlier suggestions [H 
that reconnection rates would plunge until elongated cur- 
rent layers spawned secondary magnetic islands are not 
borne out in these simulations. The rates of reconnec- 
tion approach constant values even in the absence of sec- 
ondary islands, which for anti-parallel reconnection typ- 
ically only occur transiently due to initial conditions jlOj. 
Even these transient islands can be largely eliminated by 
a suitable choice of the initial current layer width wo (a 
larger value of wo is required for the larger domains). 

A critical issue is whether the periodicity in the x direc- 
tion can influence the rate of reconnection [5[ . In each of 
the simulations we have identified the time at which the 



ion outflows from the x-line meet at the center of the mag- 
netic island. This occurs at t « 155 for the largest simu- 
lation shown in Fig.[]Ji. The plasma at the x-line can not 
be affected by the downstream conditions until t « 255, 
when a pressure perturbation can propagate back up- 
stream to the x-line at the magnetosonic speed. This is 
well after the end of the simulation. The electrons are 
ejected from the x-line at a velocity of around CAe ^> ca 
and therefore might be able to follow field lines back to 
the x-line. During the traversal time St — L x /cAe, the 
amount of reconnected flux is Vi n BoL/cA e , where vi n is 
the inflow velocity into the x-line. Using the conservation 
of the canonical momentum in the z-direction, the condi- 
tion that an electron with a velocity CAe can not cross this 
flux to access the x-line reduces to L > di(cA/vi n ) ~ 7<i, 
which is easily satisfied for the simulations in Fig. [TJ The 
fact that the reconnection rates for all of the simulation 
domains in Fig. [1] are essentially identical further sup- 
ports this conclusion. 

Also shown in Fig. [1] in the dashed lines are the rates of 
reconnection for mi/m e = 100 in (b) and rrii/m e = 400 
in (c). Consistent with simulations in smaller domains p], 
2], the rate of reconnection is insensitive to the electron 
mass. 

We now proceed to explore the structure of the electron 
current layer. Shown in Fig. [2] is a blow-up around the x- 
line of the out-of-plane electron velocity for rrii/m e = 25 
and two simulation domains, 204.8 x 102.4 in (a) and 
51.2 x 25.6 in (b), and for mi/m e = 400 in a simula- 
tion domain of 51.2 x 25.6 in (c). All of the data is 
taken in the phase where the reconnection rate and the 
lengths of the region of intense out-of-plane current are 
stationary. Reconnection forms intense current layers 
that have a well-defined length (half widths of around 
7di and independent of the size of computational domain 
for rrii/m e = 25) and then open up forming the open out- 
flow jet that characterizes Hall reconnection The 
current layer in the case of mi/m e = 400 in Fig. [2k is dis- 
tinctly shorter than the smaller mass ratio current layers 
in Fig. [2^l, b, suggesting that the length of the electron 
current layer depends on the electron mass and would be 
shorter for realistic proton-electron mass ratios. 

Shown in Fig. [3^i is a blow-up around the x-line of 
the electron outflow velocity v ex for the mi/m e = 25, 
204.8 x 102.4 run corresponding to Fig. [2^. In contrast 
with the out-of-plane current the electrons form an out- 
flow jet that extends a very large distance downstream 
from the x-line. This outflow jet continued to grow in 
length until the end of the simulation. This simulation, 
along with others at differing mass ratios, reveals that the 
peak outflow velocity is very close to the electron Alfven 
speed [8|, LUj . One might expect that because of the colli- 
mation of the outflow jet and its length, the reconnection 
rate would drop. However, this is not the case. While 
there is an intense jet in the core of the reconnection 
exhaust, the exhaust as a whole quickly begins to open 
up downstream of the current layer (J z )- The jet itself 
therefore does not act as a nozzle to limit the rate of 
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FIG. 2: (color online). Blowups around the x-line of the out- 
of-plane electron velocity for: (a) mi/m e = 25, simulation size 
204.8 x 102.4, (b) mi/m e = 25, 51.2 x 25.6, and (c) mi/m e = 
400, 51.2 x 25.6. 




FIG. 3: (color online). Blowups around the x-line for system 
size 204.8 x 102.4 with mi/m e = 25. (a) The electron outflow 
velocity v ex . (b) Momentum flux vectors, Y — p exz X- + VeyzY 
(vectors in box surrounding x-line are multiplied by 20), with 
a background color plot of | (E z + (v e x B/ c) z )/ 'E z | . 



reconnection: the rate of reconnection remains constant 
even as the length of the outflow jet varies in time. 

To understand how the electrons can form such an ex- 
tended outflow jet while the out-of-plane current layer 
remains localized, we examine the out-of-plane compo- 
nent of the fluid electron momentum equation along the 
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where v e is the electron bulk velocity, T = Pexz^ + PeyzY 
is the flux of z-directed electron momentum in the recon- 
nection plane (not including convection of momentum) 
with p e the electron pressure tensor. In Fig. [4^ we plot 
all of the terms in this equation along a cut though the 
x-line along the outflow direction from a simulation with 
rrii/m e = 100 and L x x L y = 102.4 x 51.2. The data 
has been averaged between t = 116.2 and t = 117.0. 
The electric field (black) is balanced by the sum (red) 
of the electron inertia (dashed blue), the Lorentz force 
(solid blue) and the divergence of the momentum flux 
(green). The major contributions to momentum balance 
come from the Lorentz force and the divergence of the 
momentum flux. At the x-line the electric field drive 
is balanced by the momentum transport p], [l2[. The 
surprise is that the Lorentz force, rather than simply in- 
creasing downstream from the x-line to balance the re- 
connection electric field, instead strongly overshoots the 
reconnection electric field far downstream of x-line. This 
tendency was seen in earlier simulations [12[ but there 
was no clear separation of scales because of the small size 
of these earlier simulations. Downsteam from the x-line 
the electrons are streaming much faster than the mag- 
netic field lines. Thus, in a reference frame of the moving 
electrons the z-directed electric field has reversed direc- 
tion compared with the x-line. This electric field tries to 
drive a current opposite to that at the x-line. Evidence 
for this reversed current appears downstream of the x-line 
in Fig. [2k. In spite of the strength of the effective elec- 
tric field, the reversed current carried by the electrons 
is small. As at the x-line, the momentum transfer to 
electrons in this extended outflow region is balanced by 
momentum transport. The momentum flux around the 
x-line is shown as a 2-D vector plot in Fig. [3] for the same 
run as in (a). The momentum flux has been multiplied by 
20 in the box surrounding the x-line. The data for this 
figure has been averaged between t = 172.5 and 174.5. 
The background color plot is of | (E z + (v e x~B/c) z )/E z | , 
which is > 1 where the electrons are not frozen-in. Evi- 
dent is the outward flow of momentum around the x-line 
and the much stronger outward flow of negative momen- 
tum in an extended downstream region. The momentum 
transport is so large that the out-of-plane current down- 
stream is effectively "blocked" . The force associated with 
this "blocking effect" drives the flow of the large-scale jet 
of electrons downstream of the x-line. 

We define the length A x of the inner dissipation re- 
gion as the distance from the x-line to the point where 
the Lorentz force v ex B y /c crosses the reconnection elec- 
tric field E z . At this location the effective out-of-plane 
electric field seen by the electrons reverses sign, causing 
the electron current j ez to be driven in reverse, which 
allows the separatrices to open up. Thus, the inner dis- 
sipation region defines the spatial extent of the magnetic 
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nozzle that develops during reconnection. Since the sim- 
ulations presented in this paper use artificial values of 
ra e , it is essential to understand the m e scaling of so 
that this important length can be calculated for a proton- 
electron plasma. The momentum equation of electrons 
in the outflow direction yields a steady state equation for 



-(-m eUe J = -v ez B y , 



(2) 



where v ez ~ cac- Thus, the profile of B y along the out- 
flow direction and its dependence on m e must be deter- 
mined. This profile is shown for rrii/m e = 25 (system size 
102.4 x 51.2), 100 (102.4 x 51.2) and 400 (51.2 x 25.6) in 
Fig. [4)3. Surprisingly, the profile of B y is apparently in- 
dependent of m e . Our original expectation was because 

of the continuity of the flow of magnetic flux into and 

1/2 

out of the x-line that B y ~ BoVi n /cA e ra e , where the 
outflow velocity eventually rises to cac However, since 
the electrons are not frozen into the magnetic field until 
far downstream, the expected scaling fails. To calculate 
v ex we approximate B y by a linear ramp and integrate 
Eq. ([2]). Setting the Lorentz force equation to the re- 
connection electric field, we then obtain an equation for 
A*, 



A.= [^ 

rrii 



3/8 



cE z \ 
BqCaJ 



1/2 



Bn 



doB v 



(3) 



For the three simulations shown in Fig. [4)3 the simula- 
tions yield 2.9<ii, l.Sdi and l.Odi for rrii/m e = 25, 100 
and 400, respectively, which is in reasonable accord with 
the scaling. Extrapolating to a mass-ratio of 1836, we 
predict A^ ~ 0.6^. In contrast the outer dissipation 
region can extend to 10's of di. 

We have shown that the electron current layer that 
forms during reconnection stabilizes at a finite length, 
independent of the periodicity of the simulation domain, 
and aside from transients from initial conditions remains 
largely stable to secondary island formation. Reconnec- 
tion remains fast with normalized reconnection rates of 
around 0.14. The length of the electron current layer A^ 
scales as ra^ 8 . Since the width S of the current layer 
scales with the electron skin depth c/cu pe , the aspect- 
ratio S/A x oc (me/mi) 1 / 8 . Extrapolating from our 
rrii/m e = 400 simulations to rrii/m e = 1836 should not 



significantly change the aspect-ratio and we therefore ex- 
pect the current layer to remain stable for real mass ra- 
tios. 

The structure of the current layer is important to 
the design of NASA's magnetospheric multiscale mission 
(MMS), which will be the first mission with the time reso- 
lution to measure the electron current layers that develop 
during reconnection. The length of the out-of-plane elec- 
tron current layer projects to around c/uj p i for a proton- 
electron plasma while the the outflow jet, which supports 
a strong Hall (out-of-plane) magnetic field, extends 10's 




FIG. 4: (color online). Results for simulation size 102.4 x 51.2 
with mi/m e = 25 and 100; and 51.2 x 25.6 with mi/m e = 400. 
(a) Cuts through the x-line of the contributions to Ohm's law 
for rrii/me — 100. 1 — > — m e /ev e • Vv ez , 2 — > — z • v e x B/e, 
3 -> -z • (V • P e )/(n e e), 4 -> sum of 1,2,3. (b) Cuts through 
x-line of B y for the three different rrii/m e . 
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